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by 

A.G.  Page 


SUMMARY 

Program  AcSim  was  designed  to  be  used  in  conjunction  with  a  flight 
profile  generator  program  (FPG2)  which  is  used  with  an  Intertial  Navigation 
System  simulator.  Both  AcSim  and  FPG2  were  written  in  FORTRAN-77  on  an 
ELXSI  System  6100,  AcSim  is  implemented  as  a  subroutine  of  FPG2  and  is 
called  if  the  user  selects  a  dynamic  manoeuvre  segment  to  be  performed. 
AcSim  calculates  the  dynamic  state  of  an  aircraft  for  a  sequence  of  user- 
specified  control  surface  deflections.  The  aircrafts  linear  accelerations  and 
angular  velocities  are  sent  back  to  FPG2  for  subsequent  processing. 

A  six-degree-of-freedom  model  is  used  to  simulate  the  user-defined 
aircraft.  AcSim  implements  a  State  Vector  approach  to  solve  the  non-linear 
equations  of  motion  simultaneously  using  a  4th  order  Runge-Kutta  scheme. 
An  optional  "Eigen  Analysis”  can  be  performed  which  provides  the  user  with 
longitudinal  and  lateral  characteristic  equations,  their  roots  (eigen  values), 
transfer  functions  and  unit  impulse  response  functions  for  the  given  flight 
condition. 
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NOMENCLATURE 


Symbol  Definition 

Non-dimensional  linear  accelerations,  du/dr,  dii/dr,  dwjdr 
Cl  Aircraft  lift  coefficient,  2m<j cos B„/pV2S 

y  Gravity 

/,  =  A  Moment  of  inertia  of  the  aircraft  about  Ox 

/„  -  LI  Moment,  of  inertia  of  the  aircraft  about  Oy 

I.  -  C  Moment  of  inertia  of  the  aircraft  about  Oz 

/,„  -  D  Product  of  inertia  of  the  aircraft  about  Oxy 

l^x  -  E  Product  of  inertia  of  the  aircraft  about  Ozi 

/„.  =  F  Product  of  inertia  of  the  aircraft  about  Oyz 

1.,  Non-dimensional  moment  of  inertia  about  x-axis,  A /ms2 

i„  Non-dimensional  moment  of  inertia  about  y-axis,  Il/mtj. 

i Non-dimensional  moment  of  inertia  about  z-axis,  C/ms2 
iE  Non-dimensional  product  of  inertia  about  Oxz,  E/ms'2 

h  C,J  2 

k'  -k 

L  Moment  about  x-axis  due  to  external  forces 

L,„Lr,L„,  Polling  moment  derivatives, 

it,  i ,  3h/or,ah/ar,,iL/av,f)L/3^oi./a( 

fT  Tailplanc  lever  arm  (aircraft  c.g  to  tailplane  a.c) 

/,,  Non-dimensional  rolling  moment  derivative  due  to  rate  of  roll,  I,,,  /PV„S,2 

lr  Non-dimensional  rolling  moment  derivative  due  to  rate  of  yaw,  I,r/pV„ St2 

1..  Non-dimensional  rolling  moment  derivative  due  to  sideslip,  I,,,/pV„Ss 

li  Non-dimensional  rolling  moment  derivative  due  to  ailerons,  F(/pV2Ss 

l,  Non-dimensional  rolling  moment  derivative  due  to  rudder,  Lr/pV,?Sf 

M  Moment  about  y-axis  due  to  external  forces 

Pitching  moment  derivatives, 

M„ ■„  M„  dM/dq,  aM/du,dM/dw,3M/dw,  0 M/dr) 

m  Aircraft  mass 

m, ,  Non-dimensional  pitching  moment  derivative 

to  rate  of  pitch,  M.,lpVnSt%. 

m„  Non-dimensional  pitching  moment  derivative 

due  to  velocity  increment  along  Ox,  M„/pVvS?T 
m,„  Non-dimensional  pitching  moment  derivative 

due  to  velocity  increment  along  Oz,  M,„/pV„Str 
m,„  Non-dimensional  pitching  moment  derivative 

due  to  rate  of  change  of  w,  M.i./mfr 

m, .  Non-dimensional  pitching  moment  derivative 

due  to  elevators,  M,J pV2Sfr 

rny,  rn‘H ,  ,  rn„/ii ? ,  m„,/ ip , 

m(.  ,  mj,  m,i,/in,m,l/in 

N  Moment  about  z-axis  due  to  external  forces 

/V,,,/Vr,/V„,  (foiling  yawing  derivatives. 

/Vj ,  Nr  dN/i)j>,  dN/!)r,  DN/dv,  3N/3(,  0  /V/0< 

n, .  Non-dimensional  yawing  moment  derivative  due  to  rate  of  roll,  N,,/pV, ,5s- 

nr  Non-dimensional  yawing  moment  derivative  due  to  rate  of  yaw,  Nr/pV„S.i2 

n„  Non-dimensional  yawing  moment  derivative  due  to  sideslip,  N,,/pV„S s 

n.f  Non-dimensional  yawing  moment  derivative  due  to  ailerons,  N(/pV2S* 

nf  Non-dimensional  yawing  moment  derivative  due  to  rudder,  NJpV2Ss 
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NOMENCLATURE  (cont.) 

Symbol  Definition 
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9 
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V 
Vo 

'l 

u 
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V'o 

V 
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„  XH, 
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Y„,Yr,Y„, 

Yt.Y, 

Vv 


!/r 


Angular  velocity  of  Oxyz  about  Ox  (roll  rate),  dfr/dt 
Rate  of  change  of  roll  rate,  dpjdi 
Non-dimensional  angular  velocity  about  Ox,  tijt/dr 
Non-dimensional  angular  acceleration  about  Or, 

Angular  velocity  of  Oxyz  about  Oy  (pitch  rate),  •  10/dt 
Rate  of  change  of  pitch  rate,  dq/dt 
Non-dimensional  angular  velocity  about  Oy,  dO/dr 
Non-dimensional  angular  acceleration  about  Oy,  ifO/dr- 
Angular  velocity  of  Oxyz  about  Oz  (yaw  rate),  dxli/dt 
Rate  of  change  of  yaw  rate,  dr/dt 
Non-dimensional  angular  velocity  about  Oz,  dtf/dr 
Non-dimensional  angular  acceleration  about  Oz,  d-il'/dr- 
Wing  planform  area 
Wing  semi-span 
Laplace  operator 
Time 

Time  non-dimensionalisation  factor,  m/,‘V„S 

Component  of  aircrafts  center  of  gravity  velocity  along  Ox- 

Component  of  aircrafts  c.g  velocity  along  Ox  in  the  initial  steady  state  flight 

Velocity  perturbation  along  x-axis,  U  f'„ 

Non-dimensional  velocity  perturbation  along  x-axis,  u/V„ 

Component  of  aircrafts  center  of  gravity  velocity  along  Oy 

Initial  flight  speed 

Velocity  perturbation  along  y-axis 

Non-dimensional  velocity  perturbation  along  y-axis,  v/V„ 

Component  of  aircrafts  center  of  gravity  velocity  along  Oz 

Component  of  aircrafts  c.g  velocity  along  Oz  in  the  initial  steady  state  flight 

Velocity  perturbation  along  7,-axis,  VV  tV„ 

Non-dimensional  velocity  perturbation  along  /.-axis,  w/V,, 

External  force  acting  on  aircraft  along  x-axis 
Ox  Force  component  derivatives, 

3X/3q,  dX/du.dX/dw,  3X/3n 
Non-dimensional  Ox  force  component  derivative 
due  to  rate  of  pitch,  X,,/pV„SfT 
Non-dimensional  Ox  force  component  derivative 
due  to  velocity  increment  along  Ox,  X,JpV„S 
Non-dimensional  Ox  force  component  derivative 
due  to  velocity  increment  along  Oz,  X „./pV„S 
Non-dimensional  Ox  force  component  derivative 
due  to  elevators,  X,,/pV(fS 
External  force  acting  on  aircraft  along  y-axis 
Side  force  derivatives, 

3Y/3r,  3Y/3r,  3Y/3v,  3Y/3 f ,  3Y/3< 

Non-dimensional  side  force  derivative 
due  to  rate  of  roll,  Y,,/pV„Ss 
Non-dimensional  side  force  derivative 
due  to  rate  of  yaw,  Y,/pV 
Non-dimensional  side  force  derivative 
due  to  sideslip,  Y„/pVi,S 
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NOMENCLATURE  (confc.) 

Symbol  Definition 

Non-dimensional  side  force  derivative 
due  to  ailerons,  Y,/  ,,V,‘S 
Non-dimensional  side  force  derivative 
<lue  to  rudder,  Yr/pV,;S 

Kxternal  force  acting  on  aircraft  along  z-axis 
<iz  Force  component  derivatives, 
.TZ/nq.dZ/du.dZ/dw.aZ/dri 
Non-dimensional  Oz  force  component  derivative 
due  to  rate  of  pitch,  Z^/pV„Slr 
Non-dimensional  Oz  force  component  derivative 
due  to  velocity  increment  along  Ox,  Z„/pV„$ 
Non-dimensional  Oz  force  component  derivative 
due  to  velocity  increment  along  Oz ,  Z„./p\',,S 
Non-dimensional  Oz  force  component  derivative 
due  to  elevators,  Zn/pV,?S 

Greek  Symbols 

n  Angular  displacement  of  elevator 

in  Longitudinal  relative  density  parameter,  m/,,StT 

/t,  Lateral  relative  density  parameter,  m/pSs 

t  Non-dimensional  time,  t./i 

•j,  Roll  angle 

(/■  Yaw  angle 

0  Pitch  angle 

On  Initial  pitch  altitude 

p  Air  density 

f  Angular  displacetnomf,  of  ailerons 

(  Angular  displacement  of  rudder 


1.  INTRODUCTION 

Inertial  navigation  |>lays  an  important,  role  in  air  vehic  le  navigation  due 
to  its  high  accuracy  over  short  periods  of  time,  good  performance  in  a  high 
dynamic  environment,  total  independence  of  external  systems,  and  its  immu¬ 
nity  to  jamming.  This  is  why,  even  with  the  advent  of  GRS  Navstar  (Global 
Positioning  System  Navigation  Satellite  Tracking  And  Ranging)  the  Inertial 
Navigation  System  (INS)  is  likely  to  remain  a  primary  navigation  system 
with  GPS,  or  some  other  system,  continually  calibrating  the  INS. 

An  inertial  measurement  unit  is  very  sensitive  to  errors  in  its  components 
and  their  calibration.  The  resulting  navigation  errors  grow  with  time.  Mo¬ 
ralise  different,  error  sources  (accelerometer  bias,  gyro  drift,  mis-alignment, 
etc.)  can  produce  navigation  errors  with  similar  characteristics,  INS  simula¬ 
tions  are  necessary  to  study  the  effects  of  individual  error  sources.  It  is  also 
extremely  expensive  to  test  INS’s  by  flight  trials. 

A  simulated  INS  requires  simulated  sensor  inputs,  so  the  aircraft’s  dy¬ 
namic  environment  must  be  emulated.  A  mathematical  aircraft  model  is 
required  to  obtain  a  time  history  of  the  aircraft’s  accelerations. 

A  program  exists,  FPG2  (ref.  1),  which  creates  flight  profiles  from  a 
selection  of  four  ‘ideal’  manoeuvres.  ‘Ideal’  refers  to  the  fact  that  the  air¬ 
craft’s  state  is  determined  using  simple  equilibrium  equations  rather  than 
the  flight  dynamic  equations.  As  well  as  determining  the  dynamic  environ¬ 
ment  (linear  and  angular  accelerations)  of  the  INS,  PPG 2  also  calculates  the 
nominally  true,  or  reference,  position  and  velocity  of  the  aircraft  over  the 
Parth  reference  ellipsoid. 

With  ideal  manoeuvres,  the  aircraft  proceeds  through  the  flight  by  exe¬ 
cuting  sequences  of  rotations  and  linear  accelerations,  which  FPG2  calculates 
from  the  specifications  of  each  manoeuvre.  An  ideal  manoeuvre  is  specified 
by  its  type  and  duration,  and  some  of  the  following:  heading  or  pitch  change, 
linear  acceleration  along  path,  maximum  allowed  z-axis  (See  Figure  I  for  air¬ 
craft  axes)  acceleration,  maximum  bank  angle,  maximum  angular  velocity, 
and  the  angular  acceleration  about  the  relevant  axis. 

In  order  to  improve  the  fidelity  of  the  FPG2-generated  dynamic  envi¬ 
ronment  of  the  aircraft  (and  hence  the  INS),  it  is  necessary  to  include  in 
the  simulation  a  mathematical  model  representing  the  six-degree-of-freedom 
dynamic  model  of  the  specific  aircraft.  This  document  describes  the  aircraft 
model  and  computer  program  (AcSim)  used  for  this  purpose. 

AcSim  uses  a  standard  six-degrec-of-freedom  aircraft  model  with  pilot, 
(user)  inputs  via  the  control  surfaces;  elevators,  ailerons  and  rudder.  A 
State- Vector  approach  was  adopted  to  solve  the  non-linear  equations  of  mo¬ 
tion  simultaneously  using  a  4th  order  Runge-Kutta  scheme.  This  scheme 
was  implemented  to  compliment  the  4th  order  Runge-Kutta  scheme  used 
in  FPG2.  An  optional,  so-called  ‘Figen  Analysis’,  is  also  available  which 
provides  the  user  with  the  longitudinal  and  lateral  characteristic  equations, 
eigen  values,  transfer  functions,  and  unit  impulse  response  functions.  This 
analysis  enables  the  user  to  identify  the  stable  and  unstable  modes  of  the 
aircraft  and  provides  an  indication  of  the  degree  of  stability  or  instability. 

because  FI’G2  was  written  in  FORTRAN-77  and  AcSim  was  written 
specifically  as  an  extension  to  FI’G2,  AcSim  was  also  written  in  FORTRAN- 
77.  AcSim  is  implemented  as  a  subroutine  of  FRG2  and  may  be  selected  by 
the  FRG2  user  at  any  time  during  the  simulated  (light.,  as  an  alternative  to 
any  of  the  four  ‘ideal’  manoeuvres.  If  AcSim  is  called,  the  aircraft’s  linear 
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acceleration  and  angular  velocity  response  to  the  user’s  control  surface  inputs 
are  determined  by  AcSirn  and  returned  to  F l*( J2  for  further  processing. 


2.  PROGRAM  SUMMARY 

The  ideal  manoeuvres  of  FPG2  arc  totally  independent  of  the  particular 
aircraft,  whereas  AcSim  requires  that  the  aircraft  be  fully  defined  (See  Ap¬ 
pendix  1).  If  AcSim  is  called  the  user  must  enter  the  aircraft  data,  either 
interactively  or  via  a  file. 

Once  the  aircraft  and  its  flight  condition  are  defined  the  user  is  asked 
if  an  ‘Eigen  Analysis’  is  desired.  The  eigen  analysis  consists  of  determining 
the  longitudinal  and  lateral  characteristic  equations  and  their  roots  (eigen 
values).  The  longitudinal  eigen  values  consist  of  the  Short  Period  Oscillation 
and  Long  Period  (Phugoid)  Oscillation  modes,  and  the  lateral  eigen  values 
consist  of  the  Roll,  Spiral,  and  Dutch  Poll  modes.  The  type  (real  or  complex), 
sign  (positive  or  negative),  and  magnitude  of  the  modes  give  the  user  valuable 
information  about  the  dynamic  stability  characteristics  of  the  defined  aircraft 
at  the  specified  flight  condition.  After  the  eigen  values  have  been  determined 
the  user  is  given  the  option  of  having  the  transfer  functions  and  unit  impulse 
response  functions  determined  for  each  of  the  control  surfaces. 

Even  if  the  user  doesn’t  request  the  eigen  analysis,  the  eigen  values  are 
still  calculated.  If,  and  only  if,  an  atypical  situation  occurs  (for  example,  an 
unstable  SPO  mode)  the  user  is  notified. 

After  the  eigen  analysis  is  completed  the  user  is  required  to  enter  the 
control  surface  deflection  schedule  for  the  current  sequence  of  manoeuvres. 
A  so-called  ‘manoeuvre  sequence’  consists  of  one  or  more  control  surfaces 
(elevators,  ailerons,  rudder)  being  deflected.  The  manoeuvre  is  considered 
finished  when  the  time  duration  for  required  response  output  (user-specified) 
to  the.  current  manoeuvre  control  deflections  has  elapsed.  This  data  can  be 
entered  at  any  of  three  levels,  depending  on  whether  the  user  has  saved 
manoeuvre  data  files  from  previous  runs: 

1)  Define  control  deflection  curve.  Control  deflection  curves  are  spec¬ 
ified  by  as  many  data  pairs  (magnitude  of  deflection,  time)  as  nec¬ 
essary  to  approximate  the  desired  control  deflection  curve.  'Phis 
data  may  be  saved  in  a  user-named  Uninterpolated-Manoeuvre- 
flpmipnre-File  (imMS1 2'') 

2)  Input  previously  saved  unMSE. 

Once  the  control  deflection  curve  has  been  specified,  via  I)  or 
2),  the  deflection  curve  is  assumed  linear  between  data  pairs  and 
is  interpolated  so  that  data  points  are  obtained  at  every  time 
increment,  At.  The  interpolated  data  may  be  saved  in  a  user- 
named  Manoeuvre-Sequence-Eile  (MSK). 

if)  Input  previously  saved  MSK. 

The  user  can  input  as  many  manoeuvres  as  desired  until  the  segment 
duration,  specified  in  EP(52,  has  elapsed.  Dynamic  manoeuvres  (AcSim)  is 
one  out  or  a  selection  of  five  segment  types  available  in  EI’G2. 

When  AcSim  has  the  aircraft  data  and  interpolated  control  deflection 
data  available,  the  aircraft  response  can  be  calculated.  The  linear  accelera¬ 
tions  and  angular  velocities  are  subsequently  returned  to  El’(!2  at  each  time 
increment.  The  timing  is  controlled  bv  EI>C2. 
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Figure  2  shows  the  program  summary  flowchart  of  Ac  Sim.  Figure  3  shows 
the  structure  of  subroutines  for  ArSirn.  Figure  4  is  the  group  of  subroutines 
which  performs  the  eigen  analysis. 

3.  PROGRAM  THEORY 

The  general  dynamical  equations  for  a  rigid  aircraft  in  body  axes  are 
(pi 2,  ref.  .'!)  : 
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The  aircraft  is  acted  upon  by  the  external  forces,  Y.V.Z  along  Orjhj.Or. 
respectively,  and  the  moments  of  the  external  forces,  about  Oj.Ou.n.-.  are 
/,,A/,fV  (see  Figure  I). 

Note  that  the  initial  flight  speed  vector  is  made  coincident  v/ith  Or.  so 
that  the  aircraft's  co-ordinate  system  is  defined  by  a  set.  of  stability  axes. 

Also  note  that  the  notation  used  in  this  document  and  in  the  AcSim 
source  code  for  t ho  aerodynamic  derivatives  etc.,  is  the  British  system  (used 
in  ref.  3).  A  comparison  between  this  system  ami  the  American  system  is 
given  on  p.  xxi  ref.  ft. 

3.1  Eigen  Analysis 

The  ‘Eigen  Analysis’  ultilises  the  classical  analytical  dynamic  stability 
analysis  approach.  The  above  six  non-linear  differential  equations  are  lin¬ 
earised  (on  the  assumption  that  any  disturbances  from  steady  trimmed  recti¬ 
linear  motion  are  small)  and  non-dimensionalised.  Because  the  longitudinal 
and  lateral  small  disturbance  equations  are  uncoupled  (i.e  longitudinal  and 
lateral  motion  are  independent )  two  sets  of  three  equations  are  obtained.  The 
characteristic  equations  ami  their  roots  (i.e  eigen  values)  are  determined  for 
both  the  longitudinal  and  later  .  I  sets  of  equations.  The  longitudinal  and  lat¬ 
eral  transfer  functions  and  unit,  impulse  response  functions  for  inputs  from 
each  of  the  three  control  surfaces  are  also  determined  upon  user  request.  It 
the  user  doesn't  request  the  eigen  analysis  to  be  performed  then  only  the 
eigen  values  for  the  system  are  determined  and  any  non-normal  (meaning 
not  t he  type  of  result  expected  for  a  conventional  aircraft)  conditions  are  re¬ 
ported  to  the  user.  The  theory  Tor  this  eigen  analysis  is  given  in  Appendix  2. 

3.2  Aircraft  Response  Determination 

A  State-Vector  approach  was  adopted  to  solve  the  six  non-linear  equa¬ 
tions  of  motion  simultaneously  using  ■>  fourth  order  liunge-Kutta  scheme. 

The  state  variables  were  selected  as  I',  v,  tv  and  r.  Therefore  the  State 
Vector  is: 


i  fi«'  six;  non-linear  (inferential  equations  of  motion  (l-(>)  are  re-arranged 
into  a  form  suitable  for  implementing  a  Kungo-Kutta  numerical  solving 
scheme  (Kef.  5  and  (i)  :  (Note  that  the  mass  distribution  is  assumed  constant, 
i.e  .((mays, inrrh.il) /, It  n,  and  symmetrical  about  the  roll  axis,  i.e  n) 
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where 

I,  ■<',  Hls,h,,h.lrr) 

The  forces  ( .V,  V.  '/. )  and  moments  (I..M,  \)  act  ing  on  the  aircraft  must 
be  determined  before  the  equations  of  motion  can  be  solved  for  the  aircrafts 
response. 

The  external  forces  acting  on  the  aircraft  consist  of  aerodynamic  and 
gravitational  components,  whereas  the  moments  are  only  aerodynamic  (since 
the  axis  origin  is  at  the  aircraft’s  centre  of  gravity). 

The  aerodynamic  components  are  determined  using  linear  aerodynamics 
and  the  aerodynamic  derivatives  are  assumed  constant.  Non-linear  aerody¬ 
namics  and/or  non-constant  aerodynamic  derivatives  could  be  incorporated 
if  a  special  need  arose,  but.  for  the  intended  application  both  are  unneces¬ 
sary  complications  which  would  require  a  much  more  detailed  aerodynamic 
knowledge  of  the  aircraft. 

Appendix  3  details  the  determination  of  the  forces  and  moments  acting 
on  the  aircraft. 

I’he  li  iinge-Kutta  scheme  to  solve  t  lie  equations  of  mot  ion  is  implemented 
as  detailed  in  Appendix  i. 

I'he  linear  accelerations  and  angular  velocities  are  I  hen  ret  limed  to  !•’ I ’ ( !2 
for  further  processing. 

4.  PROGRAM  VALIDATION 

Validation  of  the  AcSim  source  code  was  carried  out  by  first  validating 
the  ‘Kigen  Analysis’  section  of  the  program. 

Using  data  available  for  a  Meecbcraft  lionauza  (see  Figure  5  for  the  data 
required  by  AcSim  file  summary)  the  AcSim  Kigen  Analysis  was  performed. 
I’art  of  the  result  file  from  this  analysis  is  reproduced  in  Figure  6.  The 
characteristic  equations,  eigen  values,  transfer  functions,  and  unit  impulse 
response  functions  (for  all  .'I  control  surfaces)  were  verified  as  correct  by 
manual  calculation  of  the  results  bv  the  author  in  previous  studios. 

To  validate  the  Knnge-Kiitta  response  solution  of  the  non-linear  equa- 
t  ions  of  motion .  some  test  runs  ol  t  lie  program  were  conducted.  Kuril  control 
surface  was  deflected  separately  using  the  schedule  depicted  in  Figure  7.  To 
enable  some  comparison,  a  numerical  convolution  scheme  was  applied  to  the 
unit,  impulse  response  functions  to  obtain  the  aircraft  response  solution  to 
the  linear  small-disturbance  equations  of  motion  (see  Appendix  5).  Test 
runs  identical  to  those  performed  on  AcSim  were  conducted.  The  resulting 
aircraft  altitude  responses  for  both  solutions  are  shown  in  Figures  8a-e.  As 


l 


can  ho  soon  from  those  plots,  the  results  agree  quite  well,  with  variations 
attributable  to  differences  between  linear  and  non-linear  solutions. 

5.  ASSUMPTIONS  IN  ACSIM 

AcSim  makes  assumptions  (or  rather,  approximations)  about  the  aircraft 
itself,  which  are  : 

1.  It  has  a  longitudinal  plane  of  symmetry,  this  being  the  vertical  plane 
through  the  centre  line  of  the  fuselage.  The  mass  distribution  of  the 
aircraft,  is  assumed  symmetrical  with  respect  to  this  plane,  so  that 

lj„  0  and  /,,,  0 

Not)'  that  will  not  be  zero,  unless  Ox  and  (b  are  principal  axes  of 
inert  ia. 

1.  The  mass  distribution  of  the  aircraft  is  assumed  constant,  so  that  the 
aircraft  mass  and  inertias  are  invariant. 

.'i.  The  aircraft  is  assumed  to  be  a  rigid  body,  that  is.  any  distortion  of 
the  structure  is  ignored  (i.e  no  aeroelastic  effects). 

Linear  aerodynamics  is  assumed  by  AcSim,  and  the  aerodynamic  deriva¬ 
tives  are  assumed  constant.  This  was  done  to  avoid  complications  and  to 
minimise  the  amount  of  aerodynamic  data  required  to  model  the  aircraft. 

6.  CONCLUDING  REMARKS 

It  should  be  stressed  that  program  AcSim  is  implemented  as  a  subrou¬ 
tine  of  ri’02.  AcSim  is  an  extension  to  li'I,(;2  to  provide  the  capability,  if 
required,  of  a  more  realistic  dynamic  aircraft  response.  The  AcSim  program 
structure  and  the  implementation  of  a  dth  order  Runge-Kutta  scheme  were 
written  explicit, v  to  compliment  KI’O‘2. 

AcSim  could  be  made  into  a  stand-alone  aircraft  simulation  model,  but 
some  modification  to  the  rode  would  be  necessary.  In  particular,  the  aircraft 
altitude,  which  is  currently  passed  into  AcSim  from  l,'l><!2  as  a  direction 
cosine  matrix,  would  have  to  be  computed  in  AcSim.  Also,  some  modification 
to  the  program  structure  and  timing  would  be  required. 

The  AcSim  source  <ot)r  validation  presented  in  this  document  does  lend 
credibility  to  AcSim  results,  however  a  non-linear  solution  to  further  validate 
the  code  will  be  forthcoming.  This  will  probably  constitute  an  A.d.S.L 
(Advanced  Continuous  Simulation  Language)  model. 

AcSim  could  be  developed  in  the  future  to  incorporate  such  things  as: 
non-linear  aerodynamics,  non-constant  aerodynamic  derivatives,  or  a  ma¬ 
noeuvre  auto-pilot.  Currently,  AcSim  requires  the  user  to  fly  the  aircraft 
by  specifying  a  series  of  control  surface  deflections.  This  makes  it  difficult, 
without  many  trials,  for  the  user  to  successfully  fly  a  desired  Might,  profile. 
To  overcome  this  handling  problem,  some  kind  of  manoeuvre  specification 
rout  ine  (auto-manoeuvre  pilot)  could  be  set  up  to  provide  the  user  with  the 
option  of  flying  by  deflecting  controls  or  by  manoeuvre  specifications.  .  ,ie 
manoeuvre  specifications  would  need  to  comprise  such  items  as;  heading  or 
pilch  change,  maximum  roll  or  pitch  angle,  maximum  linear  and  angular  ac¬ 
celerations,  etc.  The  relevant  controls  would  then  be  deflected  to  obtain  the 
maximum  specified  accelerations  and  angles,  and  then  deflected  in  the  oppo¬ 
site  sense  to  reach  a  steady  condition  at  approximately  the  desired  aircraft 
stair'.  This  type  of  manoeuvre-specifying  (lying  could  easily  be  incorporated 


into  the  structure  of  AcSim  by  using  the  inano«uvrr  specifications  to  deter¬ 
mine  the  required  control  deflection  curves.  The  manual  input  of  the  control 
deflection  curves  would  then  be  by-passed.  Once  the  control  deflection  curves 
are  input  or  created,  the  rest  of  AcSim  would  run  as  usual.  Work  in  the  field 
of  manoeuvre  auto-pilots  is  currently  heing  done  at  AIM,,  and  could  be  in¬ 
corporated  into  AcSim  if  the  need  arose  to  develop  ArSiin  further.  At  the 
moment  no  future  development  of  AcSim  is  being  planned. 
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APPENDIX  1  :  Aircraft  Data  required  by  AcSirn 

The  program  requires  data  to  fully  define  the  aircraft  and  its  flight  con¬ 
dition.  The  program  will  ask  what  units  (SI  or  Imperial)  the  user  is  using, 
and  so  all  subsequent  data  must  be  consistent  with  the  tfie  specified  units. 

To  define  the  flight  condition  the  following  data  is  required  : 

Air  density,  p  (slugs/ft3  or  kg/rn3) 

Flight  speed,  V„  (ft/sec  or  m/s) 

Initial  attitude  (pitch  angle  with  respect  to  Farth),  0„  (degrees) 

To  define  the  aircraft,  the  following  data  is  required  : 

Mass  of  aircraft,  m  (lb  or  kg) 

Wing  planform  area,  (ft2  or  irr) 

Longitudinal  data  *. 

Moment  of  inertia,  or  II  (slugs. ft2  or  kg.m2) 

Tail  lever  arm  (aircraft  c.g  to  taiiplane  a.c),  lT  (ft  or  m) 

Longitudinal  aerodynamic  derivatives  (dimensional  or  non-dimensional): 

2«ci  Iff 

zu  >  zw  *  z>n  zn 

mn ,  mw ,  mw ,  rn(/ ,  rnn 

Lateral  data  : 

Moments  of  inertia  /„,  /..,  and  l,z  or  A,C,  and  K  (slugs. ft-’  or  kg.m2) 
Wing  span,  h  (ft  or  m) 

I^ateral  aerodynamic  derivatives  (dimensional  or  non-dimensional): 

,  !/,„  Mr ,  Mi ,  M< 

n„ ,  n|M  nr ,  nf 


APPENDIX  2  :  Eigen  Analysis 

A2.1  Equations  of  Motion 

The  six  non-linear  differential  equations  of  motion  are  linearised  and 
non-dimensionalised  as  given  below.  Because  the  longitudinal  and  lateral 
equations  of  motion  are  uncoupled  we  obtain  three  longitudinal  and  three 
lateral  equations.  The  equations  are  linearised  on  the  assumption  that  the 
disturbances  from  steady  trimmed  rectilinear  motion  are  small. 

Longitudinal  Equations  of  Motion  : 

l*  x"  -p  +  f  \(*\  (  **  \ 

z„  s  z„,  (I  I  j£)s  t  ~j-  tan e„  I  I  z"  I  ») 

\  /t  t  mt  a  -  n  I  m'(j  s2  —  m'f.^  /  \  ^  /  \  I1 l  m>,  / 

Lateral  E<]nations  of  Motion  : 


A2.2  Characteristic  Equations 

'I’he  program  calculates  the  characteristic  equations  (longitudinal  and 
lateral)  by  evaluating  the  determinant  of  the  matrix: 

/  s  (  A  fls  i  C  Ds  t  E  \ 

F  Qr-  I  Or  I  II  Is 2  1  .Is  I  K 

\  I,  Ms-  t  Ns  [  O  s2  I  I's  I 

'File  coefficients  (A  -Q)  are  set  at  the  bcgining  of  the  subroutines  (Fong 
and  Fat)  to  correspond  to  the  relevant  equations  of  motion. 

Longitudinal  Characteristic  Equation  : 

/S-T.„  -Z,„ 

Ai„„,,  -  del  I  z„  s  z,„  (1  (  ^  )s  I  4ft.anB„ 

\  -  /i]  -/qmts  -  /iim(„  s2  m's 

s'  1  A  As3  1  lllls7  i  CCs  l  1)1) 

Lateral  Characteristic  Equation  : 

/  *  7^s  k  (’ 

A,„,  dft  ttil„  s2  s'-’  D-s  1 

\  ^  n,,  ^-5  5^  7^  j 

A Asu  t  nils'  |  CCs'  I  IWs2  i  EEs 


A2.3  Eigen  Values 

The  roots  of  the  characteristic  equations  are  the  eigen  values  of  the  sys¬ 
tem. 

Longitudinal  Eigen  Values  : 

The  roots  of  the  longitudinal  characteristic  equation  are  then  calculated 
using  an  iterative  method  (see  Ref.  3,  pp207-9). 

For  a  ‘conventional’  aircraft  we  expect  a  solution  of  the  form: 

A . . .  (s  s,)R  -  s3)(s  sn)(s  S-i) 

where, 

.*?  i  —  rr,  \  ioj„ 

s2  n .«  ~ 

S;t  ~  (7f,  f  lUJf, 

•s*4  aV 

where 

subscript  ‘s’  denotes  Short  Period  Oscillation 
subscript  ‘p’  denotes  Phugoid  Oscillation 
V  denotes  a  complex  conjugate. 
a  -  non-dimensional  damping  factor 
w  non-dimensional  frequency 
Lateral  Eigen  Values  : 

The  roots  of  the  lateral  characteristic  equation  are  found  using  an  itera¬ 
tive  method  (see  Ref.  3,  pp21l-13). 

For  a  ‘conventional’  aircraft  we  expect  a  solution  of  the  form: 

A,„,  ,t/ts(s  -  s,)(s  •  ssHs  *s)(*  " 
where, 

s  i  <T, 

.«•>  (To 


s.\  <7;\  I  toJn 


Roll  mode 
Spiral  mode 
Dutch  Roll  mode 


A2.4  Transfer  Functions 

Transfer  functions  are  the  analytical  relationship  of  any  parameter  to  the 
input  of  any  other  input.  That  is; 

Gy  =  Transfer  Function  relating  any  parameter,  X,  to  an  input,  V'  . 
Longitudinal  Transfer  Functions  : 


A"  A"'  A" 

The  numerator  of  the  transfer  functions  are  determined  using  Cramer’s 
Rule  :  (see  Ref  4.  pp3 18-21) 
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Now,  we  have; 

Using  partial  fractions  to  obtain  a  form  that  has  a  Laplace  inverse: 
A"  (Ms3  t  (/Rs3  l  UCs  i  UP 
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The  complex  partial  fraction  coefficients  A  and  U  are  determined  by 
multiplying  and  equating  coefficients. 


Lateral  Transfer  Functions  : 


A(£<)  ,  >  A<<<>  ,.i  Al£<l 

-  -  ■  '■'If.cll'l 
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The  numerator  of  the  transfer  functions  are  determined  using  Cramer’s 
Utile  :  (see  lief  -I.  pp3 18-21) 
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Now,  we  have; 

Using  partial  fractions  to  obtain  a  form  that  has  a  Laplace  inverse: 
A;Y,,  V  As'  t VIU2  r  VCs2  h  VPs 
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The  real  coefficients  //,  and  IP  and  I)  and  the  complex  coefficient  C  are 
determined  multiplying  out  and  equating  coefficients. 


A2.5  Unit  Impulse  Response  Functions 

Unit  Impulse  Response  functions  analytically  define  the  response  of  any 
parameter  to  the  unit  impulse  input  of  any  other  input.  That  is; 

hy  =  Response  of  any  parameter,  X ,  to  a  unit  impulse  input  of  V. 

Taking  the  inverse  Laplace  transform  of  the  transfer  function  yields  the 
Unit  Impulse  Response  function  : 

where 

Ufs)  ~  £  ’If/mt  / mprWsf'l  —  1 
Longitudinal  Unit  Impulse  Response  Functions  : 

r)  =  £  'C);"""1  =  Ae •"  -h  A'e’''  t  Be’"  +  l)’e’’T 

-  2\A\e”'’  cosjw.r  +  /.A)  +  2\B\e."r’  cos(u,,r  (  /  /?) 

See  Appendix  2.1  for  the  above  manipulation  of  h*  ■ 

Lateral  Unit  Impulse  Response  Functions  : 
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See  Appendix  2.1  for  the  above  manipulation  of  A*  ,. 

(Note  :  I)  0  for  />j£t).) 

A2.6  Atypical  Eigen  Analysis  Conditions 

If  the  user  doesn’t  request  the  eigen  analysis  to  be  performed  then  only 
the  eigen  values  for  the  system  are  determined  and  any  non-normal  (meaning 
not  the  type  of  result  expected  for  a  conventional  aircraft)  conditions  are 
reported  to  the  user. 

The  atypical  conditions  include  the  following; 

-  Characteristic  equations 

-  Z‘  coefficient  in  Long-Cli  not  equal  to  zero. 

-  coefficient  in  Lat-CK  not  equal  to  zero. 

-  Real  roots  for  Short  Period  Oscillation  mode. 

-  Positive  real  part  (unstable)  of  SPO  mode. 

-  Real  roots  for  Phugoid  mode. 

-  Positive  real  part  (unstable)  of  Phugoid  mode. 

-  Complex  conjugate  root  for  Roll  and  Spiral  modes. 

-  Real  roots  for  Dutch  Roll  mode. 

-  Positive  real  part  (unstable)  of  Dutch  Roll  mode. 


APPENDIX  2.1  :  Manipulation  of  h(r) 

Mr)  -  t  /IV1" 

Aenreiujr  f-  ,u/' 

Using  :  -  cosi  |  isinz 

/i(r)  —  Aenr  (coso/r  I  i  sinu/r)  f  /l*e,Tr(tos  (--u/r)  1  i  sin  (  u >r)) 

~  enr  [  /t  (coster  f  t  sinter)  f-  /1*(coswr  tsinwr)| 

—  I  /t’Jcoswr  f  i[A  A* )  sinwr] 

,  A  A  1  lA 
now  let 

A *  A  -  iA 

so  A  \  A*  -  2 A 
and  A  A*  t2A 

therefore 

h[r)  -  2enJ  cos  ojt  A  sinter] 

-  2|ii|«rtr  r-~ -coswr  7 — rsinwr 

[\A\  \A  | 

where  \A\  is  the  magnitude  of  the  complex  number 

also,  cos  /  A  ~ 
and  sin  /  A  ■  ^|j 

[l.A  is  the  argument  of  the  complex  number  A) 
therefore 

Ji(r)  ~  2\A\eftT  [  cos  /  A  coswr  sin  /  A  sin  or  r] 

Now,  using  :  cos  x  cost/  sill  x  sin  y  ~  cos  (i  1  y) 

gives 

h{r)  -  2\A\enr  cos  (wr  I  /  A) 


APPENDIX  3  :  Aircraft  Forces  and  Moments 


The  forces  and  moments  acting  on  the  aircraft  are  represented  as  a  steady 
state  component  together  with  components  due  to  perturbations  about  the 
reference  steady  state,  taken  to  be  symmetric  wings  level  (but  not  necessarily 
horizontal)  translational  flight. 

A3.1  Forces 

The  forces  ( X,Y ,Z )  each  consist  of  an  aerodynamic  component  (X„,V„,Z„) 
and  a  gravitational  component  (X„,V,,,Z,,). 

Aerodynamic  Forces 


x„  -  X., 

„  I  uXH  I  u)Xw  I  qX.,  F  X(t) 

^  >1  ~  l  •{<■ 

I  vYv  i  pYt,  i  rYr  f  V(f) 

Zlt  ~  Zft„ 

f  uZH  I  wZv,  I  qZ.t  f  Z[t) 

'I'he  steady  state  values  are  : 

X,„, 

~  mg  sin  0<> 

Y.,„ 

-  0 

z, 

—  —mg cos  B0 

The  time  varying  components  for  this  application  are  due  to  the  control 
surfaces  only  (time  varying  excess  thrust  is  not  be  incorporated  but  could 
be  if  neccesary)  : 

x(<)  -  M')ix, 

r(t)  -  {(ok,  i  <(<m 
m  =  nW" 

Gravit.atioaa)Forc  es 

Because  the  Flight  Profile  Generator  program  keeps  track  of  the  aircraft 
attitude  via  a  direction  cosine  matrix  (DOM),  the  need  for  AcSim  to  calculate 
the  angular  displacements  of  the  aircraft  is  dispensed  with  and  the  gravity 
vector  components  are  determined  using  the  l)CM  sent  in  from  KI’G'2  at 
every  time  step. 

The  gravitational  forces  arc; 


■  X„ 

'  0 

V, 

CN 

K,I) 

n 

mg 

where  G%  is  the  transformation  (direction  cosine)  matrix  from  Navigation 
axes  (local  North, Fast, Down)  to  aircraft  body  axes.  The  direction  cosine 
matrix  from  Body-to-Nav  axes  is  calculated  in  FPG2  and  is  sent  into  AcSim, 
and  since  Off  =  (C'JJ)7  we  obtain  : 

x„-  f  n(l,3)m<7 

Y.,  ('"  (2,3)mj  -  Cjv(3,2)mg 

Z„  -  <^(3,3)mg  -  GjJ(3,3)mg 


A3. 2  Moments 

The  moments  (A, At, A)  consist  of  aerodynamic  components  At,,,  A„) 
only.  The  moments  due  to  gravity  (/„„  At,,,  A.,)  are  zero  since  the  body  axes 
origin  was  chosen  as  the  aircraft’s  center  of  gravity. 

Aerodynamic  Moments 

/..i  hvfj„  t  j>Lv  t  r/.r  t  f,(t) 

At,,  =  Af„„  t  uAt„  (  wMw  t  tt>Af„,  (  g  At.,  f  At  (£ ) 

A„  -  A„„  +  uA„  t  pA(,  I  r/tf,  r  /V(<) 

The  steady  state  values  (f,„„ ,At„„ ,A„„)  arc  zero,  since  the  aircraft  is  ini¬ 
tially  in  steady  wings-level  trimmed  (light. 

The  time  varying  components  are: 

L(t)  -  I 

M(')  - 

A(i)  C(f)A£  t  <(t)Ac 

So  the  force  and  moment  equations,  in  stability  axes,  are  : 

X  =  1  t  mgswO,,  +  (t'(t)  -  K,)X„  I  W(t)X,„  h  g(l)X.,  +  t,(().Y, 

Y  =C°(3,2)mg  1  V{t)Y„  t  r(t)Y,.  I  r(l)Yr  |  .  <(t)Vf 

Z  --  Cft(3,  I  )mg  mg  coS«„  I  (A(f)  -  V„)'/,n  I  lV(t)Z„,  I  ,(/)Z,  f  r,(t)Z„ 

L  -t  r(«)/,,.  I  r(f)/,r  i  ?(t)/,£  I  c(f)/.( 

At  =  (lt(t)  V, ,) At,,  1  u>(t)At„,  t  w[t)M,h  1  r  r,(t)  A/„ 

A  -  U(t)A„  t  ;>(t)  A,,  I  r(£) A,  I  .$(t)A£  t  .;(() A; 


APPENDIX  4  :  Runge-Kutta  Implementation 

FPG2  uses  a  fourth  order  Runge-Kutta  integration  method,  which  re¬ 
quires  the  aircraft  state  data  (linear  accelerations  and  angular  velocities)  at 
the  beginning,  twice  at  the  middle,  and  at  the  end  of  each  time  step.  In 
order  to  achieve  this  AcSirn  also  uses  a  fourth  order  Runge-Kutta  scheme  is 
implemented  as  follows  ; 

At  (  -  to 

get  A’V (,„|  using  SF„..i 
kt  -  A 


return  to  FPG2  with  .W,,.,, 

and  Ay,,.., 

At  t  to  1  ±Ai 

.sT„..i  t  k 

**1 

Rpt  AT using  .s’f'i, 

♦  4  .V)'1 

k2  At  .s'V',,  ,  iAI|i  n 
return  to  FPG2  with  .s'Vlt 

♦  I  AM,,;  aI,<^  *S  ^  it.  i  I  AM' 

At  (  to  i  t  A( 

Rot  AT,,.  AT,,..,  i  f 

k. 

AT,,.,  using  AT,,. 

♦  4  am  ’• 

ko  A/..VV  (l  ,  i  A(),.- 
return  to  FPG2  with  ,sT„ 

,  i  Arf. :•  anci  SV  v  i  An 

At  t  t„  t  Ai 

get  AT,, ...an  ,  i  k< 

get  ■'»F|,..,an  ,,s'nR  •'’'i'u.. .  .mi 

k.\  .  *  .\n 

return  to  KP<J2  with  SI 

♦  am  arifl  -s  i  if .  *  ah 

and  update  : 

1,  I,,  I  A t 

APPENDIX  5  :  Numerical  Convolution  Solution 


Once  all  the  unit  impulse  response  (UIR)  functions  and  their  derivatives 
are  calculated  the  response  to  the  given  control  deflections  (forcing  function) 
is  Tound  using  a  numerical  integration  technique:  numerical  convolution. 

The  convolution  theorem  (ref.  -I)  is: 

If  f ( t)  and  g(t)  are  the  inverse  Laplace  transforms  of  K(s)  and  O(s). 
respectively,  then  the  inverse  transform,  h(t),  of  the  product  ll(s)  K(s)(J(s) 
is  the  convolution  of  f(t)  and  g(t),  and  is  defined  by 

f 

MO  J  /(rM<  r),lr 

The  algorithm  for  numerical  convolution  is: 


c.  ,1 

■  i 

where. 

t  (i  At  lime  at  center  of  i"'  interval 
j,  value  of  reponse  at  renter  of  i"‘  interval 
f,  value  of  input,  function  at.  center  of  t"‘  interval 
h,  value  of  Milt  at  the  end  of  i"‘  interval 
i  number  of  intervals  forcing  function  is  divided  into. 

Note  that  in  the  above  algorithm  /  /(/).  and  h  h(t),  but  in  the  eigen 

analysis  theory  the  time  is  non-dimensionalised  (r  '-).  hence  we  obtain 
h  M ')■  Therefore  we  must  use  /  / ( r)  and  At  must  be  replaced  by  At  . 

We  ran  think  of  this  method  as  using  the  mils  as  influence  coefficients, 
that  is,  the  response,  say  at  t  t„,  is  the  summation  of  the  forcing  function 
at  t  !,  for  o  •  f,  <  t„,  multiplied  by  the  Hill,  with  the  mils  origin  shifted 
to  t  t,  (so  that  the  time  between  the  application  of  the  force,  at  /,.  and 
the  response,  at  to  that  force  is  /„  t,).  In  order  fot  this  to  be  done 

numerically  the  forcing  function  is  divided  into  intervals  of  width  At. 

The  elements  of  the  forcing  function  are  considered  impulses  of  magni¬ 
tude  /(t,)At  acting  at  the  renter  of  the  interval.  So.  the  response,  say  V.  at 
i  is  found  from 

.V(f„)  /(l, )*;„,(»  t,|At  .  /|t,)k;m('  Ml  At  I  *  /(U'Vu  ,(t  t.,)AC 

When  responses  are  calculated  for  parameters  with  l  ilts  derived  Irom 
other  Hills,  an  extra  term  must  be  added  to  the  results  during  the  application 
of  the  control  deflection.  To  demonstrate  this,  consider  .0 c/./r 


(i  ' 

H 

1 1) uftf'tt  ' 

1  J •'('  ')'/(<)■» 

f 

f 

| 

J  dr 

r)[f|{/).//  1 

1  'C  ('»)'/( ’  I 

FIGURE  I  :  Aircraft  axes  and  notation 
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*  Figure  from  Ref-  2,  p607 


IGURE  2  :  AcSim  Program  Summary  Flowchart 


FIGURE  3  :  AcSim  Program  Routine  Structure 


FIGURE  4  :  Eigen  Analysis  Subroutine  G 


roup 


FIGURE  5  :  Bonanza  Data  Summary  File 


***  AIRCRAFT  DATA  FILE  *** 


Air  Density  =  .002000  slugs/ft A3 

Aircraft  Mass  =  3000.000000  lbs 

Moments  of  Inertia  : 

About  Roll  (x)  axis  =  Ixx  =  1200.000000 

About  Pitch  (y)  axis  =  Iyy  =  2000.000000 

About  Yaw  (z)  axis  =  Izz  =  3000.000000 

Products  of  Inertia  : 

About  x  and  z  axes  =  Ixz  =  100.000000 

Wing  Planform  Area  =  181.000000  ft'~2 
Wing  Span  =  33.500000  ft 

Tailplane  lever  arm  =  15.200000  ft 

(aircraft  c.g  to  tailplane  aerodynamic  center) 


**  LONGITUDINAL  Non-Dimensional  Aerodynamic  Derivatives 
xu  =  -0.019  xw  =  0.056  xq  =  0 

zu  =  -0.197  zw  =  -2.12  zq  =  -0.269 

mu  =  0  mw  =  -0.306  mq  =  -0.356 

mw  dot  =  -0.0083 

x  eta  =  0 
z  eta  =  -0.4 
m  eta  =  -0.43 

**  LATERAL  Non-Dimensional  Aerodynamic  Derivatives 


yv  =  • 

-0.32 

yp  =  o 

yr  = 

0 

lv  = 

-0.06 

lp  =  -0.35 

Ir  = 

0.028 

nv  = 

0.081 

np  =  -0.039 

nr  = 

-0.123 

y  ail 

=  0 

y  rud  =  0.02 

1  ail 

=  0.02 

1  rud  =  0.005 

n  ail 

=  -0.003 

n  rud  =  -0.1 

FIGURE  6  :  Partial  Eigen  Analysis  Results  File 


***★*★★* ' 


*  ★ 

Eigen  Analysis 

Results 

*  * 

.  ** 

**- 

** 

Flight  Condition  :  Speed 

= 

290  . 

000  ft/sec 

*  ★ 

*  k 

t-hat 

S3 

889 

*  * 

*  * 

*  * 

*  * 

** 

*  A 

LONGITUDINAL  : 

*  * 

★  A 

*  * 

** 

Characteristic  Equation  = 

= 

1.0000 

sM 

k  k 

** 

+ 

7.4676 

s^3 

k  k 

A  * 

+ 

63.3017 

s*2 

k  k 

** 

+ 

1.2716 

3 

kk 

■k  • k 

+ 

1.0852 

k  k 

kk 

*  * 

★  * 

Short  Period  Oscillation 

: 

(Complex  roots) 

*  k 

(sigmaS  +/-  i*omegaS)  ** 

-  Non-dim  Damping  Factor  , sigmaS  =  -3.7247  ** 

-  Non-dim  Frequency  , omegaS  =  7.0197  ** 

Phugoid  Oscillation  :  (Complex  roots)  ’* 

(sigmaP  +/-  i*omegaP)  ** 

-  Non-dim  Damping  Factor  , sigmaP  =  -.0091  ** 

-  Non-dim  Frequency  ,omegaP  =  .  1308  ** 


LATERAL 


k  k 

Characteristic  Equation  = 

.9972 

s~5 

k  k 

k  k 

+ 

9.0336 

s'M 

k  k 

kk 

+ 

21.3877 

k  k 

kk 

+ 

94.4173 

s~2 

k  k 

kk 

+ 

1.4745 

3 

k  k 

k  k 

k  k 

k  k 

Roll  and  Spiral  Modes  :  (Real  roots) 

kk 

kk 

-  Roll  damping  factor, 

sigmal  = 

-7.8596 

k  k 

kk 

k  k 

-  Spiral  damping  factor, 

sigma2  = 

-  .  0157 

k  k 

k  k 

kk 

Dutch  Roll  Mode  :  (Complex 

roots) 

k  k 

kk 

(sigma3  +/-  i*omega3) 

k  k 

kk 

-  Non-dim  Damping  Factor, 

sigma 3  = 

-.5917 

k  k 

kk 

-  Non-dim  Frequency, 

omega3  = 

3.4138 

k  k 

k  * 

FIGURE  7  :  Control  Deflection  Schedule 


f(t)  -  control  surface  deflection  (degrees) 

/,,  peak  deflection  -  10° 
tv  —  duration  of  peak  —  Q.Ssec. 


FIGURE  8  :  Validation  Test  Runs 


Pitch.  Angle  Response  to  Elevator  T 


TIME  (soc) 


(a)  Pitch  Angle  Response  to  Elevators 


FIGURE  8  :  Validation  Test  Runs  (cont.) 


Roll  Angle  Response  to  Aileron  Test 


TIME  (sec) 


(b)  Roll  Angle  Response  to  Ailerons 


Heading  Response  to  Aileron  Test 


T 


(c)  Heading  Response  to  Ailerons 


Koll  Angle  (deg) 


FIGURE  8  :  Validation  Test  Runs  (cont.) 


Roll  Angle  Response  to  Rudder  Test 


(d)  Roll  Angle  Response  to  Rudder 


Heeding  Response  to  Rudder  Test 


(c)  Heading  Response  to  Rudder 
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